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In this article, the semi analytical solution for a fully developed time 
dependent free convective flow of a viscous incompressible fluid with heat 
source/sink in an infinite vertical coaxial cylinder saturated with porous 
material has been analyzed. The flow was induced by buoyancy forces due to 
temperature differences caused by the thermal insulation of the inner wall 
and constant heating of the outer wall. The Laplace transform technique was 
employed to transform the governing equation from time domain to 
the Laplace domain. Notwithstanding, a numerical inversing scheme known 
as Riemann-sum approximation (RSA), renowned for its precision has been 
utilized to transform the Laplace domain solution to time domain. 
The accuracy of the numerical technique employed was tested by presenting 
a comparison with the numerical values obtained using RSA, PDEPE, and 
steady state solution at large time. The effects of the various flow parameters 
on the flow formation are exhibited graphically. It is interesting to note that 
the fluid temperature and velocity increases as time passes. In addition, the 
velocity can be enhanced and minimized by gradually increasing Darcy 
number and the viscosity ratio respectively. However, the increase is seen to 
be more prominent when heat source is applied. The drag on both walls are 
seen to increase with increase in Darcy number, the reverse trend is observed 
with increase in the viscosity ratio. 
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NOMENCLATURE 

a Radius of the inner cylinder (m) 

b Radius of the outer cylinder (m) 

C p Specific heat at constant pressure (kJ/kgK) 

Da Darcy number 

g Gravitational acceleration (m/s 2 ) 

k Thermal conductivity of the fluid (W/mK) 

K Permeability of the Porous medium (m 2 ) 
Nu Dimensionless nusselt number 

Pr Prandtl number (pcp/k) 

Q 0 Dimensional heat source/sink parameter 
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R Dimensionless radial coordinate 

r' Dimensional radial coordinate 

5 Dimensionless heat source/sink parameter 

t' Dimensional time (s) 

t Dimensionless time 

T 0 Ambient temperature (K) 

T w Temperature of the hot cylinder (K) 

U Dimensionless axial velocity 

u' Axial velocity (m/s) 

V Dimensionless mass flow rate 

6 Dimensionless temperature 

GREEK LETTERS 

v Fluid kinematic viscosity 

v e ff Fluid kinematic effective viscosity 

r Skin friction 

p Density 

r] Laplace parameter 

X Radii ratio {b/a) 

y Viscosity ratio 

/? Coefficient of thermal expansion 

1. INTRODUCTION 

The study of time dependent free convective flow of viscous incompressible fluid in cylindrical 
geometry has attracted significant attention by many researchers in the past decades due to its practical 
applications in the modern industry and engineering such as in nuclear power plants, heating and cooling of 
chambers, heat exchangers, radiators, electronic cooling and air condition system. Consequently, 
the problems of free convective flow in cylindrical geometry has been solved by many researchers using both 
the analytical and numerical approach due to the complexity of the geometry and most importantly when 
the convection current is set up as a result of the internal heat source/sink in addition to the simultaneous 
heating and insulation of the walls. Several work has been devoted to examine the role of heat source/sink on 
convective flows in channel, microchannel and annulus [1-7]. 

Unsteady free convective flow of heat generating/absorbing fluid in an annulus with permeable 
walls was semi analytically studied by [1]. They concluded that that the effect of heat generation/absorption 
can be controlled by increasing/decreasing the viscosity ratio. Subsequently, reference [2] examined laminar 
unsteady fully developed free convective flow of viscous, incompressible and heat generating/absorbing fluid 
in an annulus saturated with porous material where the inner cylinder is heated constantly while the outer 
cylinder is maintained at constant temperature. They reported that in addition to the results obtained from 
their previous work [1], the heat generating fluid is desirable for optimum mass flux in the annular gap most 
importantly when the convection current is enhanced by constant heat flux. The analytical examination for 
laminar fully developed flow of viscous incompressible and electrically conducting fluid in an annulus in 
which the wall of the inner cylinder is heated or cooled either isothermally or at a constant heat flux while 
the outer cylinder is maintained at ambient temperature was considered by [3]. 

Closed forms solutions for transient fully developed free convection corresponding to four 
fundamental thermal boundary conditions in vertical concentric annulus was presented by [4-5] reported 
the analytical solution of a laminar fully developed free convection in between coaxial cylinders partially 
filled with porous material. Considerable research work has been dedicated to examine the role of heat 
source/sink on viscous flow for different physical phenomenon. In many of such problems, there may be 
plausible temperature difference between the surface and the surrounding fluid. This demands 
the consideration of temperature dependent heat source/sink which may exert strong effect on the heat 
transfer characteristics. An analytical solution for the combined effect of heat source, porosity and thermal 
radiation on mixed convection flow in a vertical annulus was presented by [6]. His findings illustrated that 
increase in the radiation parameter and heat source parameter causes an increase in the fluid temperature 
and rate of heat transfer. Recently, reference [7] analyzed the unsteady flow of a viscous, incompressible, 
electrically and thermally conducting fluid between two infinite parallel porous walls placed at y=0 and y-a. 
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It is assumed that the electrically conducting fluid is driven by mutual action of the imposed pressure 
gradient, thermal buoyancy and heat source or sink, other related problems can be found in [8-13]. 

The effect of buoyant forces on magnetohydrodynamic (MHD) free convective flow of an electrically 
conducting fluid in the presence of heat source/sink was numerically studied by [14-15] examined the effects 
of heat source/sink on MHD flow and heat transfer over a shrinking sheet with mass suction. Reference [16] 
investigated the effects of induced magnetic field and heat source/sink on fully developed laminar natural 
convective flow of a viscous incompressible and electrically conducting fluid in the presence of radial 
magnetic field by considering induced magnetic field into account and reported some interesting result, they 
obtained that increasing value of the heat source/sink parameter leads to increase in velocity. Reference [17] 
conducted a study on MHD natural convective flow of heat generating/absorbing fluid with slip effect in an 
annular porous medium. They adopted the Adomian decomposition method (ADM), a numerical scheme to 
obtain the solution of the governing equation. In the result they presented, heat generation parameter 
enhanced the temperature and velocity of the fluid in the annular gap. Moreover, slip effect parameter 
increases the velocity of the fluid. 

Employing a numerical shooting technique with a fourth-fifth order Runge-Kutta scheme, reference [18] 
carried out an investigation to study a problem of the chemical reaction and heat generation or absorption 
effects on MHD mixed convective boundary layer flow of a nanofluid through a porous medium due to an 
exponentially stretching sheet. It was found that Nusselt number is a decreasing function of the heat 
generation/absorption parameter and the chemical reaction parameter. Later on, Reference [19] examined 
the combined effects of magnetic field and heat generation/absorption on unsteady boundary-layer 
convective heat and mass transfer of a non-Newtonian nanofluid over a permeable stretching wall 
and concluded that the thermal and concentration boundary-layer thickness has higher values with 
the increasing of magnetic field and heat generation in the case of a pseudoplastic nanofluid than others. 
Other related articles by different researchers addressing some remarkable problems on the effect of heat 
generation/absorption for different fluids can be seen in [20-25]. The present paper investigates the effect of 
the heat source/sink presence on time dependent free convective flow in a vertical coaxial cylinder when 
the fluid is filled with porous material. Retardation due to Lorentz drag force opposing fluid flow is 
negligible which is therefore neglected. In addition, for effective utilization of the available heat supplied to 
the system, the inner cylinder is insulated against heat this is to retain heat within the system. 


2. MATHEMATICAL FORMULATION 

Consider a fully developed time dependent free convective flow of a viscous incompressible fluid in 
the presence of heat source/sink within a vertical coaxial cylinder of infinite length. The axis of cylinder is 
taken vertically upward in the direction opposite the gravitational force along the z'-axis, while r-axis is 
taken in the radial direction. The convective current is as a result of the constant heating at the outer cylinder 
as well as heat source/sink. The radii of inner and outer cylinder are denoted by a and b respectively as 
depicted in Figure 1. Initially, it is assumed that at time t'<0 the fluid and the two cylinders at ambient 
temperature To. At t } > 0, the outer cylinder is assumed to be heated to a temperature T w greater than that 
of the surrounding fluid and the inner cylinder which is thermally insulated. The flow is thereby setup by 
the buoyancy force due to the temperature difference of the fluid in the coaxial cylinder. Thus, under 
the usual Boussinesq approximation, the dimensionless form of the governing equations for the model under 
consideration can be written as (1-2). 



Figure 1. Schematic diagram of the problem 
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The (1,2) have been rendered dimensionless using the following non-dimensional quantities: 
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while the initial and boundary conditions in dimensionless form are given as: 

t < 0, U = 6 = 0, 1 < ft < A 
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Using U{R,r]) = f 0 °° U{R f t)e vt dt and 6(R,r]) = f 0 °° 6(R,t)e vt dt where r] is the Laplace parameter and 
r] >0. The Laplace transforms of (1) and (2) are given as: 
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The Laplace transform of the boundary conditions are given by 


(6) 

(7) 


U = 0 ; — = 0 at = 1 

dR 


U = 0; 9 =- a tR = A 


( 8 ) 

(9) 


The solutions of the Bessel ordinary differential in (6, 7) in the Laplace domain subject to 
the boundary conditions (8, 9) are given: 
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2.1. Skin frictions, Nusselt number and mass flux 

Using (11), the skin friction at the surface of the cylinders in the Laplace domain at ft = 1 and 
ft = A are given respectively as follows: 
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(S 2 -a 2 ) 

Using (10) the rate of heat transfer at the surface of the outer cylinder in the Laplace domain is given by: 

dd\ 


Nu , =-- 




= -^[^(M) - C 2 ^(A5)] 


(14) 


Thus, the solution of the mass flow rate in the Laplace domain through the annular region 
is given by: 


V = 2n j^RUdR = 2n(E 1 - E 2 ) 


(15) 


where 


E 1 = ^ (A/iOa) - ^(a)) - ^ 0*7 (Aa) - K^a)) , 

E 2 = ms L a2) ~ h ( 5 )) - C 2 {XK,{XS) - K^S))] 

It is paramount to note that the given solutions above are in Laplace domain. (10-15) is to be 
transformed to the time domain. Due to the complex nature of these solutions, the Riemann-sum 
approximation approach was adopted. Lollowing the work of [1] any function in the Laplace domain can be 
inverted to the time domain as follows: 

Z(R,t) = e ^[lz(R,r 1 )+Re(Y 1 M =1 Z(R,£ + i ™)(-ir\ (16) 

where Z represents 6, U, V or Nu as the case may be, Re refers to the real part of the summation, i = V— T, M 
is the number of terms used in the Riemann-sum approximation and 8 is the real part of 
the Bromwich contour that is used in inverting Laplace transforms. The Riemann-sum approximation for 
the Laplace inversion involves a single summation for the numerical process its accuracy depends on 
the value of 8 and the truncation error dictated by M. According to [26], the value of 8t that best satisfied 
the result is 4.7. 


2.2. Validation of the method 

To validate the numerical scheme adopted in inverting (10-15) in the time domain, we set out to 
obtain the steady state solution of (1, 2), this is achieved by computing the steady state temperature 
and velocity field analytically which is expected to coincide with transient state solution at large time. This 

can be done by taking ^ = 0 in the non-dimensional (1,2) which then reduces to the ordinary differential 
equations given by: 


, r d z u 
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R dR\ 


Da 


+ 0 = 0 


^\ + -— + S0 = 0 

dR 2 RdR 


with the given initial and boundary conditions 
U = 0; — = 0 at R = 1 


(17) 
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(19) 


U = 0; 8 = 1 alR = X 


( 20 ) 


The (17, 18) are solved under the boundary conditions (19, 20) to obtain the solution of 
the steady state temperature field, velocity field, skin friction on the surface of the cylinders, Nusselt number 
and mass flow rate. The solutions are given as: 

0(7?) = AMRO + A 2 K 0 (RO (21) 
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To further establish the accuracy of the Riemann-sum approximation approach, the PDEPE is used 
to solve (1,2) with (4, 5) as the initial and boundary condition respectively. The numerical values obtained 
from the Riemann-Sum approximation, steady state solution and the Parabolic and Elliptic Partial 
Differential Equations (PDEPE) are presented in Table 1. 


Table 1. Numerical values of the transient state velocity obtained using the Riemann-sum approximation 


method, PDEPE and that obtained from the exact solution 


Da 

T 

Riemann-sum approximation 

PDEPE 

Exact solution 

0.1 

0.2 

0.0352 

0.0353 

0.0608 


0.4 

0.0508 

0.0509 

0.0608 


0.6 

0.0571 

0.0570 

0.0608 


Steady state 

0.0608 

0.0608 

0.0608 

1.0 

0.2 

0.0541 

0.0540 

0.1128 


0.4 

0.0884 

0.0886 

0.1128 


0.6 

0.1031 

0.1032 

0.1128 


Steady state 

0.1128 

0.1128 

0.1128 


Note: (5 = 0.04, Pr = 0.71, r = 1.0, R = 1.6) 


3. RESULTS AND DISCUSSION 

In order to understand the effects of the various controlling parameters on the fluid temperature, 
fluid velocity, skin friction on the walls of the cylinder, Nusselt number and mass flow rate in the annular 
region, a MATLAB program is written to compute and generate the line graphs and numerical values for 
some selected values of the dimensionless parameter. The present parametric study has been performed over 
a reasonable range 0.001 < Da < 1.0, 0.5 < r < 1.5, — 2 < 5 < 2 where A = 2 is also taken as reference 
point. Two values of Prandtl number (Pr = 0.71) and (Pr = 7.0) were considered all through this work 
which corresponds to air and water respectively. In the course of this research work, positive values of 
S denote heat source and is specified in Figures 2-9 (a), while negatives values of S denote heat sink and is 
specified in figures. 2-9 (b ). 

Figure 2 shows the temperature profiles for different values of t across the annular gap for both 
cases of heat source and sink. It is clear from Figure 2 that the fluid temperature increases and attain steady 
state faster as time increases for both heat source and sink and drops towards the insulated wall with 
the exception of Pr = 0.71 in (a). Generally higher temperature profiles are perceived with decrease in 
Prandtl number (Pr). Velocity profiles for heat source and sink for various dimensionless parameters are 
presented in Figures 3-5. It can be seen clearly from Figures 3 and 4, that as time and Darcy number (Da) 
increases respectively, the fluid velocity is accelerated due to the fact that hot fluid particles enhance 
convective current, though the increment is significant along the heated wall and most importantly with heat 
source. 

Figure 5 presents the effect of the fluid viscosity ratio (P) on the velocity It is evident that the effect 
of viscosity ratio is to retard fluid velocity. It is seen that as the fluid viscosity. Variation of skin friction (r x ) 
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along the insulated wall are revealed in Figures 6 and 7 with respect to Darcy number (Da) and viscosity 
ratio (D) for both heat source and sink. It is found from these figures that the skin friction ( [t f) along the wall 
of the inner cylinder increases and drops as Darcy number (Da) and viscosity ratio (D) of the fluid increases 
respectively for both heat source and sink. However, the magnitude of the skin friction ( [t x ) is noticeable for 
the case of heat source. Figures 8 and 9 reveals the variation of skin friction (jf) for different values of 
Darcy number (Da) and viscosity ratio (D) along the heated wall for both heat source and sink. It is obvious 
from the figures that skin friction coefficient (rf) along the outer cylinder is enhanced with increase in Darcy 
number (Da) and decreases with increase in both Prandtl number (Pr) and viscosity ratio (D). However, 
an increase in skin friction on both walls is obtained with time. It is interesting to note that skin friction on 
both walls attain steady state faster in the case of the heat sink. It is good to also note that similar to the 
observation in Figures 6 and 7. Tables 2 and 3 reveals the numerical values for the rate of heat transfer along 
the heated wall and the quantity of fluid flowing in the annular gap respectively. It is seen that higher values 
are perceived as time is increased and heat source/sink parameter is increased that both Nusselt number 
and mass flow. 



Figure 2. Temperature profile for different values of t (5 = 2.0), 
(a) Heat source, (b) Heat sink 




R R 

(a) (b) 

Figure 3. Velocity profile for different values of t (5 = 2.0, Da = 0.1, D = 1.0), 
(a) Heat source, (b) Heat sink 
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Figure 4. Velocity profile for different values of Da (5 = 2.0, t = 0.1, r = 1.0), 
(a) Heat source, (b) Heat sink 
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Figure 5. Velocity profile for different values F (S = 2.0, t = 0.1, Da = 0.1), 
(a) Heat source, (b) Heat sink 
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Table 2. Numerical values of the transient state Nusselt Table 3. Numerical values of the transient state 
number obtained using the Riemann-sum mass flow rate obtained using the Riemann-sum 

approximation method (R = 2, A = 2) approximation method (f = 1.0, Da = 0.1,2 = 


2 . 0 ) 


T 

5 

Pr=0.71 

Pr=7.0 

t 

5 

Pr=0.71 

Pr=7.0 

0.1 

-2.0 

0.2079 

1.2078 

0.1 

-2.0 

27.3884 

28.9306 


-1.0 

0.4067 

1.2748 


-1.0 

33.0275 

29.9615 


0.0 

0.6147 

1.3422 


0.0 

38.8905 

30.9986 


1.0 

0.8328 

1.4098 


1.0 

44.9358 

32.0421 


2.0 

1.0619 

1.4778 


2.0 

51.0988 

33.0920 

0.2 

-2.0 

0.0016 

0.8365 

0.2 

-2.0 

31.0949 

31.1054 


-1.0 

0.2705 

0.9306 


-1.0 

35.8902 

32.6510 


0.0 

0.5689 

1.0256 


0.0 

39.4276 

34.2158 


1.0 

0.8985 

1.1215 


1.0 

40.5573 

35.8003 


2.0 

1.2653 

1.2183 


2.0 

37.4492 

37.4048 

0.3 

-2.0 

0.1056 

0.6554 

0.3 

-2.0 

30.3798 

33.6425 


-1.0 

0.2224 

0.7699 


-1.0 

31.7712 

35.6374 


0.0 

0.6000 

0.8861 


0.0 

27.8025 

37.6703 


1.0 

1.0422 

1.0039 


1.0 

13.3778 

39.7425 


2.0 

1.5694 

1.1234 


2.0 

-20.6518 

41.8550 


4. CONCLUSION 

A theoretical analysis on the effect of heat source/sink on unsteady free convective flow in a coaxial 
cylinder saturated with porous material was carried out. The solutions of the governing momentum 
and energy equations where derived in the Laplace domain by applying the Laplace transform technique. 
The Laplace domain solution is then transformed to time domain using a numerical inversion scheme known 
as Riemann-sum approximation (RSA). For accuracy check, the result obtained using the Riemann-sum 
approximation (RSA) approach is validated by presenting a comparison with the steady state solution at large 
time. The main findings of this work are: 

a. The fluid temperature is an increasing function of time, although the increase is more pronounced when 
heat source is applied. 

b. Velocity increases with increase in Darcy number and time. The reverse trend is observed with increase in 
the viscosity ratio. 

c. Skin frictions on both walls of the cylinder are enhanced with increase in Darcy number as time passes 
and drops as the viscosity ratio is increased. 
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